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BOUNDARY-LAYER ANALYSIS OF A PILE-UP OF WALLS OF EDGE 

DISLOCATIONS AT A LOCK 


ADRIANA GARRONI, PATRICK VAN MEURS, MARK PELETIER, AND LUCIA SCARDIA 

Abstract. In this paper we analyse the behaviour of a pile-up of vertically periodic walls of 
edge dislocations at an obstacle, represented by a locked dislocation wall. Starting from a 
continuum non-local energy modelling the interactions—at a typical length-scale of I/7—of 
the walls subjected to a constant shear stress, we derive a first-order approximation of the energy 
E-y in powers of I/7 by F-convergence, in the limit 7 co. While the zero-order term in the 
expansion, the F-limit of E^, captures the ‘bulk’ profile of the density of dislocation walls in the 
pile-up domain, the first-order term in the expansion is a ‘boundary-layer’ energy that captures 
the profile of the density in the proximity of the lock. 

This study is a first step towards a rigorous understanding of the behaviour of dislocations 
at obstacles, defects, and grain boundaries. 


1. Introduction and motivation 


Dislocations are defects in the arrangement of the atoms in a metal lattice, and their distri¬ 
bution and motion greatly affect the macroscopic behaviour of the material. Dislocations do not 
only interact with other dislocations, but also with impurities, other defects, and interfaces. Un¬ 
derstanding the behaviour of dislocations at grain boundaries and phase boundaries, in particular, 
has been the object of intensive research in academia and industry, but is still far from being 
achieved. 

The analysis of idealised pile-ups of two-dimensional point dislocations, representing the inter¬ 
sections of a system of straight and parallel dislocations with a cross-sectional plane, is a first 
attempt to shed light on this complex subject. The simplest case concerns a one-dimensional 
array of points, corresponding to dislocations in a single glide plane, forced against an obstacle 
by an external load. For this setup Eshelby, Frank and Nabarro |EFN51] pioneered a method for 
finding the equilibrium density of the dislocations, based on a polynomial representation, which is 
made possible by the specific structure of the interaction potential. This method has since been 
extended in various ways [VCMO091 IVCOAOTl IVCO07| . Other approaches to the upscaling of 
dislocations are obtained via F-convergence for phase-field energies |GM05l IFG07I ICGMllj . for 
core-radius regularized models [GLPlOl IADLGPT^ , and via homogenization for evolution equa¬ 
tions IEHIM091 iFlMl^ IMP12j . 

In this paper we study a model that lies half-way between one and two dimensions. At the 
microscopic level, the system contains a large number of periodic walls of edge dislocations with 
the same Burgers vector (see Figure |l(a)| and [RPGKO^ IGPPS131ISPPG141 IvMMPldl IHalllp . 
The coordinate system is chosen so that the walls are vertical, and represented by their horizontal 
positions a;" = (i")7o ^ [Oj oo)"+^, with ig = 0, n S N. The energy of this system is given by 



( 1 . 1 ) 


where measures the elastic properties of the medium, is the distance between consecutive 
dislocations within a wall, (t„ is an imposed shear stress, and the interaction energy potential V 


(see Figure [T(b)]) is 



1 


( 1 . 2 ) 
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The physical quantities K^, hn and tT„ are assumed to depend on n since we are interested in the 
limit behaviour of the system for n —>■ oo, and this depends strongly on the behaviour of all these 
parameters. 



Figure 1. The dislocation configuration considered in this paper. Infinite, ver¬ 
tical walls of equispaced dislocations are free to move in the horizontal direction. 
A wall is pinned at io = 0 and acts as a repellent. 


The interaction potential V is obtained by adding up the dislocation-dislocation interaction 
potentials (derived from the Peach-Kdhler force) of all the individual dislocations in a wall [HL821 
(19-75)]. Note that, although the interaction between two dislocations can be attractive or re¬ 
pulsive, depending on their relative position, the wall-wall interaction is purely repulsive. For 


this reason, the first term of the energy ( 1 . 1 ) favours configurations where dislocation walls are 


infinitely far from one another. The second term in (1.1), however, is a confinement potential, and 


models a constant shear stress acting on the system, which forces the walls against the obstacle 
at i = 0. This obstacle is a locked dislocation wall. 

As shown in [GPPS13] , the problem can be rescaled to depend only on the single dimensionless 
parameter 7 „ 


In ■■ = 


nKr, 


Roughly speaking, 7 „ is a measure of the total length of the pile-up, relative to hn- We will 
restrict our attention to the case when 1 <C 7 n ^ which corresponds to arrangements where 
dislocations are closer horizontally than vertically, while the length of the pile-up region is larger 
than the in-wall spacing /i„. 

In term of the rescaled positions Xi = 


'Ynhr 


'i—k 


the (suitably rescaled) energy ( 1 . 1 ) becomes 


g E E vMxu -»?)) + ^ E V. 

k—1j—0 z—1 

In |GPPS13] upscaled continuum models were derived from by P-convergence, in the “many- 
walls” limit n —>■ oo, for different asymptotic behaviours of the parameter 7 „. In particular, for 
the scaling regime 1 7 n ^ R, the F-limit of was proved to be the continuum energy E given 
by 

i?(/i.) = -f/ f') j p{x)^ dx + j xp{x)dx, if pL[dx) = p{x) dx, supp p C [0, oo), (1.4) 

2 7 Jg Jq 

where /r„ := 7 ^ P as n 00 . 


(1.3) 


^Our dimensionless parameter corresponds to nl 3 „ in IGPPSl^ . 
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Figurej^shows a comparison between the minimiser of the discrete energy for 7 „ = ^Jn and 
for large n, and the minimiser p* of the continuum energy E. Note that the continuous minimiser 
p* is an affine function with slope —(/^ V)~'^ 


(see Remark 


minimisers by plotting the corresponding densities; the discrete density p- 
the minimiser x" of En as 


2.2). In this figure we compare the two 
is defined in terms of 




2 /r 


i = 1 ,..., n — 1 . 


(1.5) 



Figure 2. Minimisers of and E, for n = 2^ and 7 „ = ^/n. 


The figure illustrates well the starting point of this paper. The upscaled continuum model {E) 
fits very well with the discrete model {En) in the bulk of the pile-up region; however, it fails to 
capture the distribution of dislocations at the two ends of the domain, where boundary layers 
appear. In particular, the continuum model underestimates the dislocation density at the lock, 
and understanding the optimal arrangement of dislocations at the obstacle is the fundamental goal 
of the pile-up analysis. 


Inspired by this observation, the goal of this paper is to analyse the boundary layer at the 
lock at X = 0, in the scaling regime I <C 7 n We do this by studying a F-expansion of 

the energy En |AB93[ IBC071 fBT08| in terms of the small parameter 1/^n- The zero-order term 
of the expansion is the F-limit of the energy (namely E in (1.4)), which describes correctly the 
bulk behaviour of the minimiser; the term of order 1/jn in the expansion, instead, is a first-order 
correction that captures boundary-layer effects. 

In order to capture the main features of the asymptotic expansion by F-convergence, and isolate 
them from the more technical issues related with the passage from discrete to continuum, we study 
the F-expansion of a continuum version of En- To motivate this continuum version, note that 
the discrete energies En can formally be written in terms of discrete integrals with respect to the 


measures /r„ = afj^ 


En{^ 


= \l. 


R(7„(x - y))yLn{dx)nn{dy) 


xfj-n{dx). 


Inspired by this we define the continuum energy 

nOQ pOO pOO 

■■= 7.1 I / V{l{x-y))p{dx)p{dy)+ xp{dx), 

^ Jo Jo Jo 

which now is considered to be defined for all y, € V(0,oo) - not only for sums of Diracs. This is 
the functional that we will study in this paper, in the limit 7 —> 00 . 


( 1 . 6 ) 


^This expression is formal since U( 0 ) = +oo; to make it rigorous we should remove the diagonal i = J from the 
product measure fin ® 
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1.1. Approximation of by F-expansion. The discrete energies En in (1.3) and their contin¬ 
uous counterpart Ej in (1.6) both have as F-limit the energy E in (1.41, as we show in Theorem 


2.1 


The limit energy E is therefore the first term of the approximation of E^ by T-convergence, in 
the limit 7 —)■ oo. Following |AB93) . the next term of the expansion would be the F-limit of the 
rescaled energies 

Ery(u) — imnE x 


1/7 


(1.7) 


as 7 —)■ 00 . These energies, however, do not seem appropriate for isolating the boundary-layer 
contribution at the lock a; = 0. Indeed, it is possible to show that there are two order-one 
contributions in (T^): one is due to the presence of the boundary layer (namely to the difference 
between p* and the minimiser of E^), and the other one to the approximation of E^ with E. 
The latter is actually a constant with respect to the variable p, (namely 'y{Ey{ptf) — E^p^,))) and 
does not affect the behaviour of the minimisers. Hence, to focus on the boundary-layer energy 
contribution only, we consider the functional 


') {E^{p) - E^{p^)). 


( 1 . 8 ) 


The energy in (1.8) has the advantage of measuring the difference between the minimiser of E^, 
which exhibits boundary layers, and its approximate bulk behaviour p*, in terms of the same 
energy E-^ (rather than two different energies, as in ([ill))- For this reason it is genuinely a 
‘boundary-layer energy’, since the only contribution of order one in ( 1 . 8 ) is due to the different 
behaviour of the minimiser of E^ and p* close to x = 0 . 

In order to capture the boundary-layer behaviour, we need to blow up the region next to x = 0. 
The scaling factor for the blow up arises naturally from (1.6), since the typical length-scale of 
71^(7 •) is 1 / 7 . As 7 —>■ 00 , we have 7 ^( 7 -) —t (fgV)So, which connects to the limiting energy 
E (1.4). This implies that fluctuations on the length-scale I /7 contribute to the energy difference 
( 1 . 8 ), which suggests that the boundary layer is likely to have the same length-scale. 

The boundary layer having length-scale I /7 agrees with the formal analysis of Hall ’s |Ha,ni) 
for the discrete problem with 7 „ = ^/n. Figure [^illustrates this scaling, by showing how spatially- 
rescaled densities constructed from the discrete minimisers are similar close to the lock. 



Figure 3. Rescaled discrete densities 7„p„ (see ( 1.5 )) plotted against the rescaled 
minimisers "fnx'i, for 7 „ = ^/n and for different values of n. Close to the lock the 
data points seem to lie on a single curve, confirming Hall’s prediction that the 
length-scale of the boundary layer is 0{l/^/n) = 0 (l/ 7 „). 


To capture effects on the length-scale I /7 we zoom into x = 0 by a factor 7 . Since this rescaling 
operation will be used many times we introduce a corresponding notation: given a measure p, the 
rescaled measure 7 - 1 . p is the measure 

7 _>p := 7 (x I—)■ 7 x)^p, or, in terms of Lebesgue densities, ( 7 _>p)(x) := p(x/ 7 ). (1-9) 
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This transformation zooms into the origin at rate 7 , but preserves the amplitude of the (Lebesgue 
density of) the measure. The inverse transformation is written as 7 <-/r. 

Given a measure p,, the appropriate zoomed-in version of fi describing the boundary layer is 
the measure obtained by blowing up the difference fj, — p^,, 


V :=7^[/r-p*]. 


( 1 . 10 ) 


By definition, v integrates up to zero over [0, 00 ), and it is bounded from below by — 7 _j./o*. These 
properties define the admissible class A-y in (3.2). Note that u, unlike 7 ->./r, does not see the 
mismatch of the rescaled densities outside the domain (0,1), shown in Figure]^ 

We rewrite (see Section 3.1) the energy (1.8) in terms of v as 


-I nOO poo -I poo / pOO \ 

Pii'^)-= 2 j J V{x-y)i^{dx)iy{dy)--^J \^J V{y)dyjv{dx) 

poo poo 

/ / [y {y + x))-V[y + x)\ydyv{dx), (1.11) 

Jo Jo 


1 

207 


where a = V, and all integrals (here and everywhere in the paper) denote integration over 
closed intervals. 

The main result (see Theorem |3.5[ ) is the characterisation of the limit behaviour of F-y as 7 —?> 00 . 
The main challenge in establishing this limit lies in controlling the interaction term, i.e., the first 
term in the right-hand side of (1.11). Indeed, the admissible measures v S A^ are obtained by 
blowing up probability measures, and consequently we cannot expect their total variation to stay 
bounded in the limit 7 —)■ 00 . Therefore, it is not clear how the interaction term can be defined for 
such V. For this reason, we impose non-negativity of the Fourier transform of F, which guarantees 
the interaction term to be non-negative. Moreover, this assumption allows for rewriting (see 
Lemma 3.2) the interaction term as 

1 

2./n . 


V [x - y)v{dx)v{dy) = \ [(Tvf{x)dx 


where T is the linear operator (defined in (3.13)) that satisfies T^f = V * f ior any / in a suitable 
space. We think of T as the ‘convolutional square root’ of V. This operator provides us with the 
natural framework for our main theorem: 


Theorem 1.1. The functionals Fj T-converge as j ^ 00 , with respect to the vague topology, to 
the energy F defined as 


F:A^R, F{v) ■-]- [ {Tvf{x)dx -^ / ([ V{y)dy'] v{dx), (1.12) 

2 Jr Va Jo \Jx J 

where A is the admissible class of measures defined in (3.15). In addition, a sequence with 
sup.yFy(u.y) < 00 is compact in the vague topology. 


We recall that convergence in the vague topology is defined as convergence against continuous 
functions with compact support. This topology is a natural choice, because we cannot expect the 
total variation of v.y to stay bounded in the limit 7 —> 00 , but only locally bounded. 

As a consequence of Theorem 1 1.1 1 the functional F achieves its minimum in A; since A is convex 
and F is strictly convex, this minimum is also unique. 


1.2. Approximation of the minimiser of E.y by ‘matching’. The F-convergence result in 
Theorem 0 suggests an improved approximation of the minimiser of the energy E.y at the left 
boundary of the pile-up domain, where the bulk density p* fails to describe the profile of the 
discrete density (see Figure [^. 

Denoting with u* the minimiser of F, p*(0) -I- u* is the blown-up boundary layer profile, which 
corresponds to the behaviour of the minimiser of E^ close to the lock. In view of the F-convergence 
result, we can therefore define a ‘matched’ continuous density in terms of the original, unsealed 
variables, as 

p2 •= P* + or, in terms of Lebesgue densities, pZ(.x) = p*(a;) -I- u*( 7 a;), (1-13) 









6 


ADRIANA GARRONI, PATRICK VAN MEURS, MARK PELETIER, AND LUCIA SCARDIA 


as the improved approximation of the minimiser of in (1.6). This expression appears also to 
be a good approximation of the discrete optimal density as shown in Figure]^ The agreement 
between pj and the discrete density is striking, even for a small number of dislocation walls, except 
for the free end of the pile-up region, where a second boundary layer appears, whose analysis is 
beyond the scope of this paper. 



Figure 4. Comparison between the discrete pile-up profile pn in (1.51, the ‘bulk’ 
density p* and the ‘matched’ density pJ = p* -I- 7 <-p», for n = 2 ' and 7 = i/n. 
The length of the ai-axis in the left plot equals y/aj^. 


1.3. Conclusion and comments. Theorem o gives a clear description of the boundary-layer 
behaviour of the pile-up at the lock through the minimiser of the limit energy. This theorem 
however only describes the boundary layer for the continuous energy while the ‘real’ problem 
is discrete and described by En- Still, the predictions we obtain from Theorem 1 1.1 1 are remarkably 
accurate (see Figure |4|). We now comment on some of the aspects of the result and the proof. 

Conditions on V. The F-convergence of En to E is proved in |CPPS1,^ under the assumptions 
that V is even, of class C^, integrable on K, and decreasing and convex on the positive real line. 
As a consequence of these properties, V is non-negative. A careful study of the proof in |GPPS13] 
reveals that the convexity of V can be relaxed to non-negativity of the Fourier transform of V and 
the property that V can be approximated from below in L^(]R) by continuous, positive defi nite 
function^ These properties appear again as natural requirements for the proof of Theorem 
even though the setting contains no discreteness. 

In addition, we require that V has finite first moment and that V > 0. The bound on the first 
moment is necessary to control the contribution to the energy difference in ( 1.81 which is related 
to the bulk of the full pile-up. Positivity of the Fourier transform (in addition to P > 0) is a 
technical requirement for gaining enough control on f(Tiy)^ in ( 1 . 12 ). 

Boundary layers in the discrete system. The ultimate goal of this work is to state and prove 
an analogue of Theorem o for the discrete energy En in the limit n —>■ 00 . The current results 
can be seen as a first step towards such a convergence result, in which we tackle the scaling and 
compactness issues, but side-step the discreteness. 

Unexpected ‘dip’ in the minimiser Figuregives another view on the comparison between 
the discrete and continuous densities, showing the boundary layers at both ends of the pile-up 
region in the same figure. Note that attains negative values in the region where the boundary 
layer matches with the bulk density p*. This ‘dip’ indicates that in this region the dislocation 


^To see that convexity and regularity of V implies F > 0, note that V(to) = 2 V(x) cos 2 'itujx dx = 

— V' (x) sin2'!TU}x dx', this expression can be seen to be positive by considering each period of the sine, 

and using the monotonicity of V. 
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walls tend to separate slightly farther than predicted by the bulk density. It is unclear to us why 
such a ‘dip’ occurs. 



(a) (b) 


Figure 5. Comparison between the discrete density and u* for 7 „ = y/n and 
for different values of n. Many points are omitted to enhance readability. See 
Section |4^ for a detailed description of these plots. 


In conclusion, Theorem 1 1.1 1 together with the numerical illustrations predict accurately the lo¬ 
cal aggregation of dislocation walls at a lock. This aggregation is not recovered by the existing 
continuum models for dislocation pile-ups |EBG04[ IHallll IGPPS13] in the corresponding param¬ 
eter regime 1 ^ 7 „ ^ n, because these models only predict the affine density prohle of the bulk 
behaviour. Hence, these continuum models break down at a mesoscopic length-scale (given by 
0(1 /jn) for 7n ^ ^/n, based on numerical computations) near the lock, which is larger than 
the microscopic length-scale of neighboring dislocation walls, but smaller than the macroscopic 
length-scale of the full pile-up. As a result, these continuum models underestimate the stress 
around the lock. 

Theorem 1 1.1 1 provides a way to fix this error by taking the effect on the mesoscale into account. 
Our results imply that the length-scale of the dislocation density is too coarse to describe the 
behaviour of dislocations close to grain boundaries. However, by taking into accounts effects on 
an intermediate mesoscopic length-scale, we can obtain an accurate description of the dislocations 
distribution in the pile-up region. We believe that such a mesoscopic description can be used to 
improve continuum models for the dislocation density at grain boundaries. 


The plan of the paper is as follows. In Section]^ we prove the convergence of to E (i.e., 
the zero-order approximation of Ey); we then prove the main result. Theorem in Section 
In Section we discuss the agreement between discrete minimisers and the matched continuum 
density p2 (1.13) by means of numerical computations. Finally, Appendix [A| and are devoted 
to some technical steps in the proofs of results in previous sections. 


1.4. Assumptions and notation. Although the expression (1.21 is the inspiration for this paper, 
and although we use (1.2) in the numerical calculations, the convergence results hold for a broader 
class of functions V. In Sections and we make the following assumptions on V: 


(VI) V : K —>• K is even, and decreasing in (0, oo). 

(V2) V G L^(R) has finite first moment, i.e., |V(x)| dx, /jj |ai| V(a:) da; < oo. 

(V3) V has positive Fourier transform, which satisfies 


V9 G if2’~(k). 


(V4) There exists a sequence (14) C C'b(M) of functions with non-negative Fourier transform 
such that 14 t ^ pointwise a.e. and \\V — 14|jii(R) —>■ 0 as fc —)■ oo. 
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Note that by (V1)-(V3), V is non-negative and |lt^||oo = ^(0) = fgV(x)dx < oo. We need 
the technical condition (V3) to define the ‘convolutional’ square-root of V by the operator T as 
introduced in Theorem 0 (see Section [3.1[ ), and we rely on (V4) to rewrite the interaction term 
in (1.111 in terms of T. 


Here we list some symbols and abbreviations that we use throughout the paper. 


a 

A, A-y 

En 

E^ 

E 

I E{f) 

E-Hf) 

E, 

E 

h-y 

Al([0,oo)) 
, v~ 

P([0,oo)) 

P* 

P* 

T 

X 

II • IItv 

For fi G M{[0 


= /o°° ^ = I /h ^ in the case of IQ])) 


’ dx 


admissible sets for F and F-y 
transforms of p by scaling space by 7 > 0 
discrete energy 
continuous energy at finite 7 
(zero-order) limiting energy 
Fourier transform of /; T{f){uj) = /(w) := Jjj 
inverse Fourier transform of /; 
first-order energy at finite 7 
limiting first-order energy 
auxiliary function in the expression of F-y 
one-dimensional Hausdorff measure 
signed Borel measures on R with support in [0,00) 
positive and negative part of a measure v € Af ([0, cio)); > 0 
non-negative Borel measures on [0, 00) of mass 1 
minimiser of E 

rescaled version of p*; p* := 7-»p* 

‘convolutional square root’ of V; f = V * f 

Hilbert space on which T : X ^ X is bounded, and A <Z X' 

L'^-norm, for 1 < g < cx) 

total variation norm of a signed Borel measure. 

00 )), we define f := /j^ / dpi. 


(3.151, (3.21 


(1.91 


(1.31 


(1.61 


(1.41 


(Till 


( 1.121 


(3.51 


(2.11 
Sec. 4.1 

and (2.21 

(3.1f 

(3.12 

1, ( 

1 

3T3f 


2. Zero-order Gamma-convergence 

In this section we derive the zero-order term of the F-expansion of in ( |1.6[ ) in powers of 
1 / 7 . This term is the F-limit of as 7 —>■ 00 , in the narrow or weak topology (i.e., convergence 
against continuous and bounded functions). Not surprisingly, it turns out to be the same as the 
F-limit of the discrete energy En in ( 111 , namely the continuum energy E defined in ( |1.4[ ). This 
is proved in the following theorem. 

Theorem 2.1. The T-limit of E^ os 7 — > 00 with respect to the narrow topology is given by 

noo poo 

I / p{x)'^ dx + / xp{x)dx, if p{dx) = p{x) dx, 

Jo Jo 

otherwise, 

where a := V. In addition, sequences with bounded E^ are compact in the narrow topology. 

Proof. Compactness, i.e., tightness, is a consequence of the inequality E.y{p) > xp{dx). The 
proof of the liminf inequality follows from a simple modification of the proof of the discrete-to- 
continuum convergence in |GPPS13l Theorem 7]. 

As for the limsup inequality, we first reduce by density arguments to p G Cc([0,oo)). Then 
we take as recovery sequence the constant sequence p.y := p. Note that Vj(x,y) := ^V{'y{x — y)) 
converges vaguely to 2a'H^cD, where D is the diagonal of the first quadrant, 

D := {{x,y) G (0,oo)^ |x = y}. 

Then, we conclude the limsup inequality from 


E{p) = \ JO 
-\-oo 


lim / 

7-*-oo7o 


Vj (x, y)p{x)p{y) dx dy = 2a 


/ p{x)p(y)l-L^\-D{dxdy) = 2a / p 

/o ^0 


□ 
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Remark 2.2. We recall that the minimiser of E is given by 


'’•<*> I - 1 ; 


/X 4 
E{p*) = gVa, 


hence the rescaled version of p* is, according to the definition (1.9), 




3. First-order Gamma-convergence 


( 2 . 1 ) 


( 2 . 2 ) 


In this section we prove Theorem o (see Theorem |3.5[ ), which states the T-convergence of 
in (1.111 to F in (1.12). 

3.1. Preliminary results. Rewriting the energy. Here we show that under the blow-up trans¬ 
formation ( 1 . 10 ), which transforms p into v, the energy E.y converts into the expression for Fly 
in (ll lll). 


Setting jF := we write 


= A 


r >00 fOO 


/o ^0 


V{x — y)fP(dx)fP(dy) H— / xjp{dx) =: E^{fP). 

7 Jo 


V{x - y)n{dx)i^{dy) + J V{x - y)p^{y/j)dy + i^{dx). (3.1) 


Note that fp = 7 . By defining as in (1.10) v \= jp — we write, by using Fubini’s theorem, 
Fy(u) := 7 {E^{p) - E^{p^)) = E^ifl'^) - Ejip*) 

^ pOO nC 

Jo 

For the sake of notation, we omit the superscript 7 . The measures v obtained via the rescal¬ 
ing ( 1 . 10 ) belong to the admissible class defined as 

.4y, := S A 1 ([ 0 , 00 )) : J |i/| < 27 , J u = 0 , iy{dx) >—pi,{x)dx^ . (3.2) 

Whenever convenient, we extend u to M by zero without changing notation. Note that by (3.2) 
supp iy~ C supp = [0, 2'yy/a] for each v G A-y. 


and ( 2 . 2 ), 


(3.3) 

Next we rewrite F.y{v). By using the explicit expression of p* in (2.1) and the fact that V is 
f 1 

y{x-y)p^{y/-f)dy- 


even, the integrand of the second term in (3.1) can be cast into 

/»cm 

X 


>0 


7 


I PlVa. 


Plpa 2 ^ r 

/ y{x-y)ydy +-— / V{x-y)xdy 

JO ^*^7 Je 

1 I r I r 

= -r y{x-y)dy +-— {x - y)V{x - y)dy + -— / V{x-y)ydy 

Jo 2a7 Jr 207 JH\(0,27y/o) 

1 1 1 f °° 

=y{x-y)dy +-— / y(.x-y)ydy +-— / V(x - y){y - 2-fyPi)dy 

V ® Jo ^®7 J-00 ^®7 J2-iy/a 

1 /■°° 

= 2 Ja-p/ V{y)dy + h^{x), 

V ^ J X 


= —j= / V{x-y)dy- — 

pa Jo 




(3.4) 


where we have set 


1 1 

h^{x) := -^V * ((a;A0) V (t - 27 Ja)) = — J [y {y + {2'y^/d-x)) -V{y + x)]ydy. 


(3.5) 
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Substituting (3.4) in (O and using the fact that u = 0 , we obtain 

^ poo poo poo 

F^{v) = - / {V * i'){x)iy{dx) — / g{x)v{dx) + / hj(x)i'(dx), (3-6) 

2 Jo Jo Jo 

1 r°° 

where g{x) := ^ / V{y)dy. (3.7) 

yO' Jx 

Note that, due to assumptions (VI) and (V2), g e C'b([0,oo)) n L^(0,oo) is non-negative and 
non-increasing. 

Properties of the auxiliary function . In the following lemma we derive some useful estimates 
for h^. The qualitative behaviour of is illustrated in Figure]^ 


7 defined in 


is non-decreasing and satisfies the bounds 


\hP < C- 
7 

in 

[ 0 , 27 v^] , 

(3.8) 

IGI < 

in 

[V7, 2jy/d - v^] , 

(3.9) 


for some 17 ( 7 ) —>■ 0 as 7 —>■ 00 . 


hy{x) 



Proof. The function h.y is non-decreasing, since it is the convolution of a non-negative function 
with a non-decreasing function. 

To prove (3.8) we note that by the monotonicity of hj, for all x G [0, 2^^/a]^ 

1 11/’°° 

— C— < h.y{Q) < h~f(x) < h.y{2"f\/a) < C-, with C = — / V(y)ydy. (3.10) 

7 1 2 a Jo 

As for ( |3.9[ ), we estimate 

1 f°° 1 f°° 

h^{2lVa-y/^) =< -— / V {y P y/^)y dy < -— / V{z)zdz, 

2^7 Jo 207 

which implies (3.9) by assumption (V2). □ 


The ‘convolutional square root’ ofV. In |GPPS1.^ Appendix] it was shown that, for the special 
case of ( 1 . 2 ), 


V{lu) = / V(s)e-"""‘" dx = 


1 


2ui sinh( 7 r^a;) 


cosh(7r^a;) — 




sinh(7r2w) 


We observe that V is even, strictly positive on K, and p := \Pv £ IF^’°°(M) is strictly positive. 
However, since V decays to zero at infinity at rate l/|a;|, the function p{uj) decays at rate 
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and fails to be in L^(IR). Therefore, the inverse Fourier transform of p may not be a function IT, 
and it is not clear whether there exists a function IT : M —>■ K such that V = W * W. 

In this paper we also want to consider any V satisfying conditions (V1)-(V4). We therefore 
circumvent these difhculties by defining ‘convolution with IT’ as an operator. We leave the details 
of the construction of this operator to Appendix and focus here on its key properties. We 
consider the operator 

T:L^(E)^L^(E), Tf:=J^-^(pf). (3.11) 

Since p^ = V, we obtain 

T^f = T-^{p^f) = T-^{Vf) = T * /, 

which motivates us to think of T as the ‘convolutional square root’ of V . However, in view of the 
expression of F in Theorem we need to extend the domain of T beyond L^(IR). It turns out 
that a natural space for this extension is the dual of the Hilbert space 

:= {/ e H\E) : x^f{x) S L^{E)}, (/, ^)x := J /V' + J /<^(1 + x^). (3.12) 

The space X is chosen to be as large as possible to require little regularity on p, but small enough 
such that X' contains the domain of F. Lemma 3.2 and Lemma 3.2 below quantify this 
statement. 

Furthermore, Lemma [3.2| lists the properties of T and X which we need in our proof of Theorem 
O In particular. Lemma |3.2[ ( p^ applies to G Ay for all 7 > 0 (by choosing / = p*). 
Consequently, the first term of Fy{vy) in (3.6) equals the first term of F{vy) in (1.12). The proof 


W := 


in Appendix [A| 


of Lemma 3.2 (^iiib is a corollary of Lemma 3.4 We prove the remaining statements of Lemma 3.2 


Lemma 3.2. The linear operator T defined by (3.11) and the space X defined by (3.12) satisfy 
the following properties: 

(i) T is a symmetric bounded linear operator from L^(IR) to itself, from X to itself, and from 
X' to itself, where we define 

T-.X'^X' by {T^,f)~{f,Tf); (3.13) 

(ii) For all ^ G X' it holds that = V * 

(Hi) j u € A4([0, 00 )) : v{dx) > —p,t{0)dx, sup u“'"([a;, a; + 1]) < 00 1 C X'; 

(iv) For all v G At([0,oo)) for which there exists f G H^(E) with xf{x) G L^(E) such that 
V + f > 0 is a finite measure, it holds that 


f {Tv)^{x)dx = f {V * v){x)v{dx). 

JR JO 


Domain A of the limit energy F. We now give the complete definition of the limit functional F: 

F:A^E, F{v) = f I {TvY(x)dx — / g{x)v{dx), (3-14) 


= \ [ {Tvf{x)dx- [ g{x)v{dx), 

^ JR JO 


where g is given in (3.7) and where the admissible class A is 


A'.= Iv G A1([0, 00 )) : v[dx) > —p*(0)da;. 


supu''‘([a;,a; + 1 ]) < 00 , 

ai>0 


Tv G 


(3.15) 


As for Ay, we extend u S A to K by zero whenever convenient. By ‘Tu G L^(E)’ we mean that the 
linear operator Tv : (X, || • II 2 ) —>■ M is bounded; Lemma 3.2 (i]),(pil|) imply that this requirement is 
meaningful because of the lower and upper bounds on v in A. 

While for finite 7 the set Ay consists of measures with uniformly bounded total variation, in 
the limiting set A this bound has vanished, and a priori it is unclear why a sequence (u'*') C Ay 
with bounded energy Fy(v'^) could not acquire infinite total variation in the limit. Instead, our 
compactness result in Section 3.3 implies v G A, which gives a local bound on \v\, but no global 
bound. Maybe minimising sequences have bounded total variation, but we have not investigated 
this. 
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3.2. Key estimates. Lemma [373] states two estimates which are central to the proof of the com¬ 
pactness statement in Theorem o Lemma |3.4| shows how the local bounds ot v G A provide 
control over the linear term f gv. 

Lemma 3.3. Let rj > 0. For any bounded interval [r, s] C [0,cx)) with s — r > rj there exist 
constants C > 0 and Cjj > 0 (decreasing in g) such that for all j > 0 and v G A^^ the following 
estimates involving \ •=Tu hold 


< CriiVs - r||A ||2 + s-r), 

/ O rO 

TX= V*v<C{\\X \\2 + l). 

-oo j —oo 


(3.16) 

(3.17) 


Proof. It is not restrictive to prove the claim only for u with Tv G L^(]R), the bounds (3.161 and 
(3.17) being trivial otherwise. 


Since V is positive, we obtain 


y * = 


V{x — y)v~^{dy)dx > 


r J 0 


V{x-y)dx] v+{dy) > p+([r, s]), 


where is a strictly positive constant increasing in g and converging to 0 as ry 0. Since we 
assume that Tv G L^(]R), it follows from Lemma 3.2 (^ that T^v G L^(IR) and ||r^p ||2 < C||Tp|j 2 . 
We then get the desired upper bound (3.16) from Lemma 3.2 1© by 

Cr^v^ifr^s]) < f V * v^ = f V*v+ f V * v~ 

Jr Jr Jr 

pS PS POO 

= / T'^vF / / V{x — y)v~ (dy) dx 

Jr J r Jo 

< y/s-r ||T^p ||2 + (s - r)p*(0)||y||i 

< CV^WTvh + {s- r)p*(0)||y||i. 

For proving the second estimate (3.17), we use ( |3.16 ) with s — r = g = 1. Since V is positive, 
even and increasing in (—cx), 0), we have the estimate 


/•O poo 


V * V < 


< 


' —oo J 0 


V{x — y)v^{dy)dx 
V{x — k)v~^{[k, k + l])dx 

k=o 

^ pOO 

<C{\\Tvh + l ) J 2 

k=0''^ 

where the sum in the right-hand side is finite due to assumption (V2). □ 

Lemma 3.4. There exists a constant C > 0 such that for all v G At ([0, oo)) it holds that 

gv < Csup |p|([a;, X-I-1]), where g is as in (3.7), and (3.18) 

x>0 

) 

fv < C'||/||ASup|u|([x,x-b 1]) for all f G X. (3.19) 


x>0 


Proof. Let ip G {g} U X. We set Ik ■= {k, k + 1) for A: G N, and estimate 
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We estimate the right-hand side of (3.20) iov ip = g and p ^ X separately. Since g G L^(0,oo) is 
bounded and non-increasing, we obtain 

oo pQo 

XI ll5lU“(/fc) < 5 ( 0 ) + / g<oo, 
k=0 


hence (3.181. Next, to prove (3.19), we set (p = f £ X, and fix fc S N. By continuity of /, there 
exists Xk S Ik at which p{xk) = P- Therefore, by the Fundamental Theorem of Calculus, we 
obtain 

f{y) = fpk) + r{f{x)ydx = [ /' + 2 r //' for all y G 4. 

Jxk J Ik IXk 

Then, by using the Cauchy-Schwarz and Young inequalities, we estimate, for fc > 1, 

ll/lli“(/fe) ^ /^ + 2||/||L2(/^^)||/'||i2(j,^) < -j^\\Pf\\p(Ii^) + J^\\Pf\\L'^(Ik)\\f'\\L^iIk) 


'Ik 


LHIk) 


II/' 


\L^{Ik) 


'' ¥' 


(3.21) 


Note that f G X implies {bk)k C We obtain the desired estimate, by using (3.21) and the 

Cauchy-Schwarz inequality on sequences, from 


El 

/c=l 


0 ^EtV^- 






k 


OO 1 

'Ll 


.fc = l 


A)2 


S'-. 

.fc=l J 


<C\\f\\x, 


and from ||/||l~(/o) < CH/Hx, by the first inequality in (3.21). □ 

3.3. Main result: T— convergence. The following theorem is the main result in this paper. 

Theorem 3.5. (First-order T-convergence). Any sequence {v-)) C with bounded energies 
Fjpj) is pre-compact in the vague topology. Furthermore, Fly T-converges to F os 7 —)■ oo with 
respect to the vague topology, i.e.. 


for all v.y ^ v vaguely, 


liminf Fly(j 2 y) > Fiv), 


7 —>-oo 


for all 11 G A there exists ^ v vaguely such that limsupFly(^ 7 ) ^ Pp)- 


(3.22a) 

(3.22b) 


7 —>-oo 


Proof. Compactness. Let (vj) C yly be such that Fy(zzy) < C. We split the proof into several 
steps. 

Step 1. Fp-f ) is bounded. We claim that jy hyi/y is bounded from below. This bound implies, 
together with Lemma 3.2 10 applied to (jZy), that 


poo 

F(yP = Fjpj) — / hyi2y < C. 

Jo 


(3.23) 


From estimate (3.23) we conclude in particular that finiteness of F-^ pp implies finiteness of F{uj). 

We now prove the claim. Since hy is non-decreasing (Lemma 3.1) and by (3.3) supp C [ 0 , 27 y 7 ], 
we deduce the bound 

poo poo p2^y/a poo poo 

/ /lyZ2y = / hyp — / /lyZ2y > /ly (0) / ^^y ” hj (2^ ^/a) / 

Jo Jo Jo Jo Jo 

where the last inequality follows from (3.2) and (3.10). 

Step 2. Convergence of v-f and T. From Step 1 we have 

1 

C > F{i^P = PXPI-gv^, 


(3.24) 
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where := We rewrite the second integral in the right-hand side of (3.24) in terms of A-,, 
by 

poo ^ poo / poo \ ^ poo / poo \ 

J \J ^ y yj Viz + x)dzj v^{dx) 

= ^ [ ( r v(z- x)dz) lyj(dx) = ^ r = ^ r TX^, 

y/a Jr VJ-oo / VO y/a 

where exchanging the order of integration is justified by V G (K) and || P 71 | tv < 7- In conclusion, 

from (3.24) and Lemma 3.3 we obtain 

C > F(v,) = 2 IIA^II^ - ^ TX, > illA^II^ - CdlA^lb -f 1); 

hence (A^) is bounded in L^(]R). Consequently, by (3.16), v^{[r,s\) is uniformly bounded in 7 for 
any bounded interval [r, s] C [0, 00 ), with a bound that only depends on s — r. Therefore, along a 
subsequence. 


V vaguely, and X~^ = Tv^ ^ X in L'^. 


(3.25) 


Step 3. Characterisation of the limits v and X. By the basic properties of the vague topology, 
V e A4([0, 00 )) satisfies i'{dx) > —p»{Q)dx and sup 2 ,>o v+([a;, x-I- 1]) < 00 . Consequently, by 
([iiit, V S X'. We claim that converges weakly-* in X' to v. If so, the continuity of 


Lemma 


3.2 


T (Lemma 


3T^([I|) and (3.25) imply that A = Tv S L^(ffi), and thus v & A. 

To prove this claim, we take f G X arbitrary and a large M > 0 fixed and write 

poo pM poo 


= fv^+ /V 7 . 


/o 


10 


IM 


The first term in the right-hand side converges to /v as 7 —>■ 00 because of the vague con¬ 
vergence of Vj to V. We now show that the second term is infinitesimal for large M. Indeed, by 
Lemma 13.41 and Lemma 13.31 


= 


IM 


{fX(M,oo))v^ < C'||/X(M,oo)||ASUp|v7|([x,X-f 1]) 

tc>0 


= C'll/X(M,oo)ll V sup(v| -k )([x, X + I]) 


a:>0 


< C||/x(m.oo)I|a(I 1A7||2 + 1 + p*(0)), (3.26) 

where X(m,oo) denotes the characteristic function of the interval (M, 00 ). Clearly the term in ( 3.26| ) 
is infinitesimal for large M, as claimed. 


Lower bound. To prove (3.22a), we observe from (3.23) that it is sufficient to show that 

poo 

liminf / hjV^>0, and liminf F(v.^) > F(v). (3.27) 

7—>-00 Jq 7—>-oo 

We start with the first inequality in (3.27). First of all, if we write h^ = hfj — h~ and , 

we have 




h+v^ - 




the other terms being non-negative. Therefore the first inequality in (3.27) follows if we prove 
that 


pOO poo 

limsup / h1fv~ = 0 , and limsup / h~vff = 0 . 

7—)-oo Jo 7—foo Jo 


(3.28) 


We note that, since suppv.^ C [ 0 , 27 y 7 | and v (dx) < pif(fS)dx, 


poo p2^y/a pz^y/a 

/ < p^{Q) / 

Jo Jo Jo 


277/5 
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The bounds (3.8) and (3.9) on entail 




h^= h+ 
Jo 


■27V“-VT 




h+ 


/ hi:; < C'ct(7 ), 

\J1 


(3.29) 


(3.9) and Lemma 3.3 we have 


which converges to zero as 7 —>■ oo. Analogously, since (supph..^) n [0, 00 ) C [0, 27 -^ 0 ), by (3.8) 

• 27^0 


= 




h~vt + / Kvt 


7 7 


hZvt 


C 


' vA 
, cr(7) 


< — / + C ■ 

1 Jo J Jvz 

C 

< —— + C(t(j) —)■ 0 as 7 
V7 




J2^y/a-^t 

(J /■27%Ai 

H- / •'-f 

1 J2^y/a-y/^ 


V, 


,+ 


00 . 


(3.30) 


The bounds (3.29) and (3.30|) above imply (3.28), and hence the first inequality in (3.27). 


We now prove the second inequality in (3.27). From (3.14), (3.25) and Step 3 of the compactness 
proof, we obtain immediately that 

liminf ||ru-y ||2 > \\Tv\\l. 

7 —>-oo 

Therefore, to conclude the proof of the liminf inequality, it is sufficient to show that 


gu^ 


gv as 7 —)■ 00 . 


Since g is continuous, it follows from the vague convergence of Vry to v that gv^ —>■ gv as 

7 —> 0 for any M > 0 fixed. The remaining part of the integral is arbitrarily small for M large 


enough, uniform in 7 , by Lemma 3.3 and Lemma 3.4 


Upper bound. Let v £ A. We first show in Steps 1 and 2 that by the usual density argument 
it is sufficient to prove (3.22b) for v e L^(0,oo) with bounded support and v > —p,(0). 


Step 1. Approximation of v in L^(0, 00 ). We leave the proof to Appendix [ b| 

Step 2. Truncation and retraction. For v G ACi L^(0,oo), we set 

Vn ■■= (1 - ^)uL[0,n]. 

It is easy to see that S A n L^(0, 00 ) satisfies the desired bound > —p*(0), and that 
converges both in the vague topology and in L^(0, 00 ) to v. Con verge nce of F(i/„) to F{v) as 
n —)■ 00 follows from L^-continuity of F : L^(0,oo) —>■ K (by Lemma 3.2 (^ and g G L^(0,oo)). 

Next we construct a recovery sequence {v.y) for v G AO L^(0,oo) with support in [0,M] and 
V > —p*(0). We define 

rM 


Vj{x) := v{x) — a^pt{x), 


1 

CT.y := — / V. 

7 Jo 
-2 


It is clear that v.y converges to v in both the vague and L^-topology. We conclude that v.y G A^, 
since v.y =Q and v.y > —p* for sufficiently large 7 . 

It remains to prove that F.y{v.y) -G F{v) as 7 —>■ 00 . We use (3.23) to write 


plVl p 

= F{iy^) + / h^iA — CTj 

Jo Jo 


2^/07 

-ylJ (T-y / . 


(3.31) 


By using again continuity of F, we obtain that F{vj) -G F{y) as 7 —>■ 00 . With the bound on 


the estimates 


given by (3.8) we show that the remaining two terms in the right-hand side of (3.31) are small by 

pM /-7 n2y/a'f 


h.yV 


< \JM—\\v\\2, 

7 


and 




< ( 7 ^,(21/07)—p*(0). 
7 
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□ 


4. Numerics 


In this section we present some numerical results which illustrate the scaling properties of the 
boundary layer and the degree of fit or misht between the predictions of the continuum limit 
energies E and F, and the minimisers of the discrete energy i5„. The numerical simulations 
presented in this section are all obtained for the specihc choice of the potential given in ( 1 . 2 |. 


4.1. Rescaling of the density. The density v appearing in the boundary layer energy in (3.6 1 
is defined by blowing up the pile-up domain by a factor 7 . This behaviour is confirmed by the T- 
convergence result. We now illustrate some formal calculations and numerical results that suggest 
a similar scaling for the discrete densities. 

As predicted by the asymptotics in [Halil) . Figureillustrates that the size of the boundary 
layer region at the lock is 0{l/^/n) in the case 7 „ = ^/n. In Figure we have computed the 
dislocation wall positions cc” by using Newton’s method to solve \/En{x) = 0. To define a discrete 
density in terms of the rescaled coordinates 7 ^ 2 ;”, we mimic the transformation ( 7 n)-> (defined 
for the continuum densities in (1.9)): 


Pn{x) ■■= Pnixljn), X G { 7 „a 


i = 1, 


■-1}- 


(4.1) 


Although Figure [^illustrates the size of the boundary-layer prohle, the slopes of the densities 
at the right-end of the domain are different. This difference is a result of the n-independent 
bulk behaviour of the unsealed densities which converge to the bulk prohle p* (cf. [GPPS131 
Thm. 7]). 

For this reason it is more instructive to compare the densities := Pn — p*, which we do in 
Figure The good agreement between the prohles of for different values of n and different 



. n = 2® 

□ n — 2® 

X n = 2*^ 

□ 

■ 

□ 




_^ , °^° ° ° oB Q q 1 , 


(a) 


0.7797 - X 


X In = 
n 7n = s/n 
" ln=rdl^ 




51°. .q, „ in,. 


(b) 


Figure 7. Plots of the rescaled density for (a) 7„ = and various values 
of n; (b) n = 2 ® and various scalings of 7 „. 


scaling of 7 „ (conform 1 <C 7 n ^ r) suggest that also in the discrete case the density must be 
rescaled as in (14.1|). 


4.2. Rescaling of the energy. From the analysis in this paper one could try to infer that also 
for the discrete energies i?„ in (1.3) the right scaling to analyse the boundary layer at the lock is 
l/ 7 „, namely that \En{x^) — En{p^)\ oc 1/jn- Here, p" are the discrete dislocation wall positions 
generated by p*, i.e., 


P* = Vi '■= 2^0 (^1 “ \/l “ “ 


rVi 


solves 


p* — 


n 
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We investigate this question by means of numerical computations, for 7 „ being a power of n. 
Setting an ■= \En{x'2) — En{p'^)\, and assuming a power-law relation also for namely n~P, 
we try to determine p by computing 


Pn ■ = 


log an - log a 2 n 
log 2 


for n = 


n3 n4 

z , z ,.. . , z , 


(4.2) 


for different expressions of 7 „. The results are shown in Table The data in the first two columns 
support our conjecture that « log„ 7 „ for sufficiently large n, although the third column seems 
to contradict it. This disagreement could be due to a much slower convergence to the correct 
value, or to the presence of a second boundary layer at the free end of the pile-up (see Figure [^, 
whose energy contribution interferes with the contribution of the boundary layer at the lock. At 
this stage we can only guess as to the reasons for this discrepancy, and we will return to this issue 
in a future publication. 


n 


log„ In 


1 

4 

1 

2 

3 

4 

23 

0.550 

0.901 

0.931 

2^ 

0.384 

0.757 

0.833 

2^ 

0.299 

0.642 

0.745 

2® 

0.263 

0.567 

0.675 

2^ 

0.250 

0.525 

0.625 

2® 

0.247 

0.505 

0.590 

29 

0.247 

0.496 

0.567 

210 

0.247 

0.494 

0.552 

2“ 

0.248 

0.493 



Table 1. Exponents p in (4.2) for a power law fit of \En{x^) — A„(p")| oc 
for different values of n and 7 „ = The lower-right value is missing 

because of memory limitations in MATLAB. 


4.3. Minimiser of the continuum energy E. In this section we analyse the minimiser of F, 
namely we consider the minimisation problem: 

minF(u), 


where F and A are defined in (3.14) and (3.15). We assume in addition that the minimiser is 
in L^(0, oo) and satisfies the strict inequalitjQ inf u* > —/9*(0), and derive a first-order necessary 
condition for minimality for u*. We consider variations u* -I- t/r S A for p € Al([0,oo)) with 
bounded support and bounded Lebesgue density, and |t| sufficiently small. Then 

0 = —F(u* + tM)|t=o = y {V*u^)p-J gp, 


and since p is arbitrary, this implies that 

1 


V *v^= g = —j=V * X(_oo.o] 

V Cl 


a.e. on (0,oo). 


(4.3) 


To approximate the solution to (4.3) numerically, we follow the method in Section 10.5 of 

(4.4) 


[HeaOlj . We approximate by a step function of the form 

AT 

u*« y] Xixu > 


^Note that the numerically calculated minimisers (Figure are consistent with this assumption. 
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where Xi are the unknowns, while N and the sets li are chosen beforehand. Figure [^suggests to 
take small intervals li = {ai-i,ai) close to 0, and larger intervals in the region x > 1. We choose 
the intervals li so that \Ii\ = CV, and then tune C, b and N such that oi = 3 • 10“®, 200 S /jv, 
and the interval Ij containing the element 1 has length \Ij\ = 0.1. This results in C = 2.727• 10“®, 
b = 1.1, and N = 141. 

After having fixed and N, we substitute the approximation (4.41 in the Euler-Lagrange 
equation (4.3), and we solve for Xt the linear system 


XI * XiXVj) = 

i—1 ^ 


).o])(2/j) for yj = 


®7-l 




(4.5) 


This amounts to an asymmetric Galerkin (collocation) approximation, in which the solution space 
is approximated by piecewise constant functions, and the test space by sums of delta functions 
supported on the midpoint of each interval. 

Note that we can identify {V * X/J(yj) = fo; and that we have explicitly 


V = Li2(e"2^) - yln(l - e"^^). 


where Li 2 is the polylogarithm, 


OO t. 






The solution to the linear system (4.5) results in the profile of u* (4.4) as depicted in Figure]^ 



Figure 8. 
the x-axis. 


Numerical approximation of u* (4.4). The intervals A are depicted on 
The circles indicate the coordinates (i/i, u»(j/i)). 


In Figure we compare u* to the discrete density constructed from the min¬ 

imiser pn of En- In order to make the figure cleaner, instead of plotting the (mathematically 
correct) piecewise-constant version of Figure]^ we plot the affine interpolation of the coordinates 
(j/D u*(j/i)). The figure shows a very good agreement between and u* in the boundary layer 
region. As we have seen from Figure [7 the profile of Vn depends neither on n nor on the scaling 
of 7 „. Therefore, together with Figure ^ we conclude that u* is the boundary-layer profile of the 
dislocation density at the lock x = 0. 

On the other hand, the agreement between and u* is not very good outside of the boundary 
layer region, as shown in Figure]^ where we plot the densities on (l,2^nVa) on a logarithmic 
scale (remember that (0, 2"fn-\/0') is the whole of the pile-up domain, after blowing up the variables 
by 7 „). In Sectionj^we proved that F-converges to F with respect to vague convergence, hence 
minimisers of the energies Fly converge to u*, minimiser of the limit energy F, only on bounded 
domains. It is therefore not surprising that the densities u„, the discrete analogue of the minimiser 
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Figure 9. Comparison between the discrete density for n = 2^^ and 7„ = ^/n 
and the minimiser u* of F. Only the interval (0,1) is shown, which corresponds 
to the region (0, l/ 7 „) for the unsealed densities. 


of Fly, do not agree with u* on the whole pile-up domain, which is of order 7 „ in the blown-up 
variables, and hence growing in size with n. 

From Figure we can actually infer that there is another ‘continuous’ boundary layer at the 
right-end of the pile-up domain, since there is a large number of particles deviating from the bulk 
distribution p* (i.e., with density deviating from 0). We expect the profile of this boundary 
layer to resemble the bulk behaviour of the optimal discrete density of the energy (1.31 for 7 „ = n. 
We refer to |GPPS13l Theorem 8] and to Figure 6 in the same paper for the analysis of that case. 
The analysis of the boundary layer at the free end of the pile-up region is however beyond the 
scope of this paper. 


Finally, as noted in the introduction, (or u*) curiously attains small negative values before 
matching the line y = 0. As this ‘dip’ is present in the numerical approximations of both and 
u*, it does not appear to be a numerical artifact. We remain with questions regarding how this 
‘dip’ relates to the choice of V. Does every convex decreasing V with finite first moment have 
such a ‘dip’ ? Are there choices of V for which u* is decreasing? Does V have more than one ‘dip’ 
(the numerical solutions for and u* are too coarse to give a speculative answer), and how does 
the answer depend on the choice of VI The answers to these questions are beyond the scope of 
this paper. 
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Appendix A. Properties of the operator T 


Here we prove Lemma 3.2 in more generality through Lemma |A.2| The generality comes for 
free after introducing the Hilbert spaces in (A.l I, which are naturally related to the Fourier 
transform. 
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First, we recall in Lemma |A.1| several basic properties of the Fourier transform and of the 
complex Hilbert spaces defined by 


:= {/ 


= <f&H 

k 


‘^(K;C) : x^f{x) G L2(M;C)}' 
fW(l)Wf f {x)(j){x) dx 


for all k, j G N. 


e=o • 


In case the functions are real-valued, we set 


Xk,j := {/ 




: x^ f{x) G 


(A.l) 


In particular, we note that X defined in (3.12) equals Xi^ 2 , and that L^(M) = D 5(M) 

for all k,j G N, where 5(M) is the space of Schwarz functions. The analogous statement for 
complex-valued functions holds as well. 

Lemma A.l (Basic properties of T and Xkj{C)). 

(i) F{x !-)■ x^ f^^\x)) = i-A- for all k,£ G N; 

(ii) If fJ, G V{M.), then Xfi is a positive definite, bounded, uniformly continuous function; 

(Hi) IF is a symmetric, invertible, bounded linear operator between Xkj{C) and A^- fc(C) for all 
k,j G N; 

(iv) For k G N+, Xk,kiC) and X^^k are Banach algebras with respect to both the pointwise product 
and convolution; 

(v) Let k G N+, / G Xk,k{C) and £, G Xk,k(Cy. Then F{f * 0 = /C Xk,k{Cy. 


Lemma A.2. The linear operator T defined by (3.11) satisfies the following properties. 

(i) T maps L^(M) to itself; 

(ii) For all k G N and j = 0,1,2, T is a symmetric bounded linear operator that maps Xkj to 
itself, and pi := {ip i—)■ pp) is a symmetric bounded linear operator that maps Xj^k to itself; 

(Hi) For all k G N and j = 0,1,2, the operator pi is symmetric, bounded and linear, and maps 
X( J, to itself. Furthermore, the extension of T to X'f. j given by 

for^G Xy^, f G Xk,j, (T^J) := {^,Tf) = {F-\pO,f) 
is well-defined as a symmetric bounded linear operator mapping X'f, ^ to itself; 

(iv) For all ^ G A 2 2 , it holds that = V * and ||T ^||2 = J n® elements of [0,oo]; 

(v) For all v G Al([0,oo)) for which there exists f G Xij such that n + f > 0 is a finite measure, 
it holds that 


[ {Tv)^{x)dx = f {V * v){x)v{dx). 
Jb Jo 


C L° 


), it holds that pf G L^(IR;C), 


Proof of Lemma [ A. ^ We now prove the claims (i)-(v). 

(i) . Let / G L^(M) be arbitrary. Since p G 
and thus Tf G L^(IR;C). Moreover, since f{—oj) = f{uj) and p is real-valued and even, we also 
have that {pf){—uj) = (p/)(a;). Hence, Tf = F~^(j)f) is real-valued. 

(ii) . Since T = F~^(pI)F, Lemma l^diiit implies that it suffices to prove that pi := (p i-G- pp) 

is a symmetric bounded linear operator mapping Xj^k to itself. Symmetry is obvious. We prove 
boundedness by showing the existence of some C > 0 such that ||p/||j,fc < for all / G Xj^k- 

We obtain such a constant C from the following computations, which rely on p belonging to 

for 0 < £ < j G {0,1, 2} we have 


moreover 


iPff^W < e (0 

^ i=0 ^ J ^ i=0 

II^V||2< |b||oo||a.Vl|2< Wooll/IU 
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(Hi). It follows directly from the previous step that pi and T are symmetric bounded linear 
operators on X( ^ and ^ respectively. Since p is an even function, we obtain the representation 
formula for T( with ^ from 

inj) = {^Tf) = {i,X-^{pf)) = ilpf) = ilpf) = (pII) = {T-^{p^j), for all / £ 


(iv). By (V 2 ) we have that V £ L^(IR) C i. Then, Lemma A.l implies that T^f = V f 
for all / £ Xk^k for any k £ N-|-. By using the extension of T in the previous step, which requires 
fc < 2, we obtain ^ for all ^ £ X2 2- 

Next we prove that as elements of [ 0 , 00]. To avoid confusion, we introduce 

a different notation for the extension of the L^-norm to linear functionals. For any 77 £ X2 2, we 
define 

||i7||l2(k)' := sup w^e[0,+oo], 

feX2,2 II/II2 


which is consistent with the usual definition of || 
over all / £ L^(IR)) because ^2^2 is dense in L^( 
The assertion follows from step (Hi) by 


|i2(Ry (in which the supremum is taken instead 
i). Likewise, we interpret as ||p^|| 


L2(R)/- 


IIK|Il2(R)/ = II J'T^||l2(h)/ = ||T^||i2(H),, 

where the Fourier transform does not change the value of the L^(IR)'-norm. 

(v). It is not clear how /(F * u)u is defined without using the asserted structure that / £ Xi^i 
and /i:=u + /> 0 isa hnite measure. Therefore, we define 

y'(F*z2)u = I (V^{p-f)){p-f):= J{V^p)p -2 y'(F*/V + J {V * f)f. (A.2) 

We claim that the right-hand side is well-defined in (—00,00]. Indeed, since / £ Xi^i, we have 
that V * f = T^f £ Xi^i C C't,(M). Moreover, since /i is a hnite measure, the second and third 
integrals are hnite. Since V and /i are non-negative, we hnd that F * /r > 0 , and thus the hrst 
integral is dehned as an element of [ 0 ,oo]. Moreover, since ||F * < ||F||i||/r||Tu < 00, Fubini 


applies to J(V * f)p, = f (V * p)f, which shows that the dehnition in (A. 2 ) is consistent with the 
case in which f(V* p)p < 00 . 

As in step (iv), we interpret | 1 Fu||| 2 (r)/ £ [0,oo]. Here we rewrite it in terms of 

/ and /i. Whenever convenient, we regard /i as a bounded linear operator on C';i(K), dehned by 
(/r,/) := Jg^fdp. Since C{,(]R) D Xi^i, it follows that p £ Cf,(K)' C X[^^, and thus Tp £ X[ i is 
well-dehned by step (Hi). 

From Tv = T/r+ (T/, • )2 we obtain that Tv £ L^(IR)' <;=^ Tp £ L^(M)'. If Tv £ L^(IR)', then 
we set /jy and as the Riesz representatives in L^(IR) of Tu and Tp respectively. We obtain 

rHli 2 (Ry = IIMI^ = 11 / 7211 ^-2(/^,r/)2 + ||T/||2 = \\Tp\\l 2 (^y- 2 {Tp,Tf) + J (F*/)/. (A.3) 


If Tv ^ L^(IR)', then we claim that the right-hand side in (A. 31 equals -l-oo. Indeed, Tv ^ L^(IR)' 
implies Tp ^ L^(K)' and thus \\Tp\\L2(^^y = -|-oo. It remains to check that \{Tp,Tf)\ < 00. Since 
p £ X[ i and / £ Xi^i, we have that Tp £ X[ i and Tf G Xi_i, and thus {Tp,Tf) is indeed hnite. 
In conclusion, (A. 3 ) holds for any v as in the assumption. 


It remains to prove that the right-hand sides of (A. 2 ) and (A. 3 ) are equal. The equality of the 
second terms in both expressions follows from 


{Tp,Tf) = {p,T^f)= J{V^f)p, 
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and thus it remain to show that f {V * = ||r^||^ 2 (Ry in [0, oo]. We prove this statement by 

establishing the following results: 


y (y * < oo 


\\T^i\\ 




^mII 


LH 


< oo 




(A.4a) 
(A.4b) 


are 


We first prove (A.4a). Given any n G N, let Vn be as in (V4). Since Vn,V G L°^ 
non-negative, they define (semi-)inner products on L^(IR) by 

(/, V)v := J {V * f)(p = J Vf0, and (/, ip)v^ := j (K, * f)ip = j Vnf0. 

We construct Hilbert spaces as the closure of L^(R) with respect to the corresponding norms || • ||v 
and II • ||y„: 

Hv := and Hv„ := L2(M) 

Since it is given that f (V ^ < oo, we obtain /r G Hy- Because Vn < V and ^ > 0, it holds 

that / {Vn * < J{V * and thus fj, G Hv„ ■ Therefore, we obtain 

/ {Vn * = sup 2 / {Vn * - j {Vn* ■ 

J _ J J 

Next, we fix (p G Xi^i, and use 14 —>■ H in L^(IR) and 14 t ^ pointwise together with the Monotone 
Convergence Theorem to pass to the limit n —>■ oo in the previous expression, to get 

j {V * = lim / (14 * > lim 2 / (14 * ^p)fJ. - {Vn * 

J n—>-oo J n—^oo J J 

= 2 J{V*^)^l- J{V*ip)ip = 2(T/i, T^) - J {Tipf. 

Taking the supremum over p G Xi^i, we obtain 

[{V*fi)fj.> sup 2{Tp,(j))- f (jV , 

J (/)GTXi,i L j 

where TXi.i := {Tf : / £ It follows from /? > 0 that TXi^i is dense in Xi^i. Therefore, 

J{v*fi)p> sup 2{Tn,p)-Jp^ =\\Tfi\\l 2 (^^y. 

Next, we prove (A.4b). Since it is given that L^(IR)' 9 Tp = X~^{pll), it follows from Ji G 


C'b(M;C) (by Lemma A.l (^) that pp is a measurable function, and thus we can write 

J \pp\^ = I V\p\^. 


\Tf^\\ 


L2(R)' 


(A.5) 


We introduce two regularisations: we regularise V again by 14, and we set := Gg * p G 
L^{R) n C'f,(M), where is the scaled Gaussian satisfying 

G^{x) := exp(—Tiex^) ^ 1 vaguely as e —>■ 0, and Gg = e~^l'^Gii^ ^ dg as e —)■ 0. 

First, we use 14 € Gf,(M) in combination with ® ^ ^ (g) ^ as e —>■ 0 to deduce that 

{Vn* = // Vn{x - y) d{pg p^){x, y) // Vn{x - y) d{p p){x, y) = {Vn*p)p. 

J J JK2 j Jts ,2 j 

Then, we use that Pe > 0 and 14 < H to estimate 

J {Vn * Pe)lJ-e < J {V * Pe)Pe = J 4|/i^P = J V Gl\p\^ . 

Since f 1 pointwise on K, we obtain from the Monotone Convergence Theorem and (A.5) that 

f VGl\p\^^ [ V\p\^ = \\Tp\\l 2 ^^y. 
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Summarizing our arguments, we have 


J (\4 * n)n = Im J (14 * < \\Tn\\l2(j^y. 


We conclude by taking the supremum over all n. 


□ 


Remark A. 3 (Assumptions on V needed for Lemma A.2). For T to be a bounded operator from 
X to itself (cf. Lemma |A.2 


ii^), more regularity than p £ L° 


is required. In the proof we 


showed that p £ 


The proof of Lemma A.2| 


is sufficient. 

is technical, and we need (V4) to complete it. Instead of the 


requirement p > 0, it is sufficient to have \{p = 0}| = 0. However, we need the requirement p > 0 
elsewhere (for constructing a recovery sequence) to control the tail of u. 


Appendix B. Technical step in the proof of (3.22b) 


i.e., we show that it is sufficient to prove (3.22b I 


Here we prove Step 1 of the proof of (3.22b) 
for u G A n L^(0, oo). 

First we show by density arguments that we can restrict ourselves to those v £ A with support 
away from 0. For any u £ A, we consider Vc ■= for any c > 0, where t^{x) := x + c acts as 
a horizontal shift. It is clear that Uc converges vaguely to u as c —>■ 0, and that Uc satisfies the 


upper and lower bounds in the definition of A. Since by (3.11) T commutes with t^, it follows 
that ||ri/c ||2 = ||Tu ||2 < 00 . Hence, Vc £ A, and to complete the density argument it remains to 
show that limsup 2 _,,Q J gvc = J 9^- By splitting the integrals as follows 

/ p pM poo 

gu - j gvc= j [gix) - g{x + c)] u + / [g{x) - g{x + c)] u, 

J Jo J M 

it is clear that for any M > 0 the first term in the right-hand side converges to zero, because g 


is continuous, and the second term is small in M by Lemma 3.4 and Lemma |3.3| by a similar 
argument as in Step 3 in the proof of Theorem |3.5| 

Next, we take v £ A such that suppu C [c, oo) for some c > 0. We consider the approximation 
4 := 'ipg * u, where is the rescaling of a carefully chosen mollifier ip (constructed below), and 
its restriction := 4i-[0,oo). We choose ip as the Fourier transform of 

(1 - a|cLi|)^(3a|a;| + l), if |a;| < 1/a, 

0 otherwise. 


ip(uj) := 


where a > 0 is chosen such that 4 = 1. The function ip is constructed in |Wen041 Table 9.1] to 
have compact support and C^-regularity while being non-negative and positive definite |Wen041 
Thm. 9.13]. Then, by Bochner’s Theorem, we also have that ip > 0 is positive definite. By basic 
properties of the Fourier transform we find that 0<ip<l, 0<ip<l, '0(1) = '0(1) = 1, and 
ip,ip £ AI3 3 C X. For the scaling in e, we set 


^■,(u = 4(j). 


and thus 0e(a;) = ipieio). 


We note that ipg <5o narrowly and 0£ t 1 pointwise as e —>■ 0. 

To check that = UeL[0, 00 ) is admissible, i.e., € A, we observe from ipe > 0 having unit 

mass that the bounds on lyf are satisfied. We show that Tv^ £ L^(IR) later in this proof. 

Now we prove the vague convergence of to u. Let M > 0, and let tp £ (7b([0,oo)) be an 
arbitrary function with support in [0,M]. Since ip^ is even, we rewrite 

p]\4 pA/[ pAd /•Af+1 poo 

/ V’Ue = / = / (^(0e*u)= / {ip^: * ip)v + / {lps*p)v. (B.l) 

Jo Jo Jo Jo Jm +1 

Since ipe * p converges to p uniformly in M as e —^ 0, the first integral in the right-hand side of 

since 


(B.l) converges to J ipu. The second term in the right-hand side is bounded by Lemma 
Ipe * <p £ X, and moreover small in e since p is supported in [0, M]. 


3.4 
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We observe that the second term of the energy F{y^) converges to the corresponding term of 
the energy -F(u), i.e., 


gve 


gv as e —)■ 0, 


by Lemma 3.4 and g G (^^([OjOo)). It remains to show that G L^(R) and that the first term 
of th e energy F(Pe) satisfies the bound limsup£_j.Q ||Tj/e ||2 < ||^rj/|| 2 . Since V > 0, from Lemma 

implies that 


A.2 


we deduce that 9 e Li^^(U;C). Then, since ijje G ^ 2 , 2 , Lemma 


A.l 


FVg = l/'e * u = ^^9. 

We conclude from ijjs having bounded support and 9 (]R;C) that G L^(M;C), and thus 

Finally, we show that limsupg_,,Q ||Tp £||2 < ||Ti/|| 2 . To this aim, we split 

||TPe||2< ||r9,||2 + ||T(9,L(-c3o,0))||2. (B.2) 

that 


We obtain from i/’e < 1 and Lemma 


A.2 


||T9,||2= / I-|9,| = / !/>, I-|9p< / y|9p = ||Tp||2. 

»/ ]R IR IR 


It remains to show that the second term in the right-hand side of (B.2) is small in e. Since 


ij}" G L^{R), G C'b(M), and therefore ipe{x) < Cex Since suppu C [c,-l-oo) with c > 0 

and sup 3 ,>Q v{[x,x + 1]) < 00 , it follows that for x < 0 


\9e{x)\ ^ 


(x — zy 


'{dz) < C'e^ 


fc=0 


(—X -|- fc -I- c)^ 


Hence, ||9e||i2(_oo_Q) —0, and thus the second term in (B.2) is small in £. 
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